Belle Preprint 201 1-9 
KEK Preprint 2011-7 



Observation of 5 — )■ pAD'^ at Belle 
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We study meson decays to pAD^*'^ final states using a sample of 657 x 10*' BB events collected at 
the T(45) resonance with the Belle detector at the KEKB asymmetric-energy e+e^ collider. The observed 
branching fraction for B^ — > 'pAD^ is ( 1 .43^g^5 ±0.18) x 10^^ with a significance of 8. 1 standard deviations, 
where the uncertainties are statistical and systematic, respectively. Most of the signal events have the pA mass 
peaking near threshold. No significant signal is observed for B^ — >■ pAD*" and the corresponding upper limit 
on the branching fraction is 4.8 x 10^^ at the 90% confidence level. 

PACS numbers: 13.25.Hw, 13.30 Eg, 14.40.Nd 



Since the first observations of baryonic decays of B mesons 
by ARGUS ll and CLEO H], many three-body baryonic B 
decays have been found Although the general pattern of 
these decays can be understood intuitively from heavy b quark 
decays 101, many specific details cannot be explained by this 
simple picture. 

Using a generalized factorization approach, Ref. Js]] 
predicts rather large branching fractions (^10^^) for the 
Cabibbo-suppressed processes B pAD^*\ The branch- 
ing fractions of other related baryonic decays such as B'^ ^• 
ppD° HQ], B° ^ ppK*° Mi, B- ^ ppK*- H m and 
B — >■ ppn^ istl are used as inputs in such estimates because 
baryon form factors entering the decay amplitudes are diffi- 
cult to calculate from first principles. The expected values 
of the branching fractions for B pAD^ and B — > pAD*" 
are already within reach with the data sample accumulated at 
Belle. 

Nearly all baryonic B decays into three- and four-body fi- 
nal states possess a common feature: baryon-antibaryon in- 
variant masses that peak near threshold. This threshold en- 
hancement is found both in charmed and charmless cases fl- 
A similar effect has been observed in J/\\t ^ p'pj decays 
by BESJpl [il and CLEO I^L but is not seen in ^ 
p'pTi^ lUlll and Y(15) — > ppj II14I1 . One of the possible ex- 
planations of this phenomenon suggested in the literature is a 
final state A^A^ interaction ifTsll . 

In this paper, we present results on the B 'pAO''*^^ de- 
cays in order to test the factorization hypothesis and study the 
pA threshold enhancement effect. 

The data sample used in the study corresponds to an inte- 
grated luminosity of 605 fb^', containing 657 x 10*' BB pairs, 
collected at the T(45') resonance with the Belle detector at 



the KEKB asymmetric -energy e+e^ (3.5 GeV and 8 GeV) 
collider lfl6ll . The Belle detector iITtIi is a large-solid-angle 
magnetic spectrometer that consists of a silicon vertex detec- 
tor (S VD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel-like ar- 
rangement of time-of-flight scintillation counters (TOE), and 
an electromagnetic calorimeter (ECL) composed of CsI(Tl) 
crystals located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. 

The selection criteria for the final state charged particles 
in B pAZ)" and B pAD*" are based on information 
obtained from the tracking system (SVD and CDC) and the 
hadron identification system (CDC, ACC, and TOE). The pri- 
mary and D*^ daughter charged tracks are required to have a 
point of closest approach to the interaction point (IP) that is 
within ±0.3 cm in the transverse (x-y) plane, and within ±3.0 
cm in the z direction, where the +z axis is opposite to the 
positron beam direction. Eor each track, the likelihood val- 
ues Lp, Lk, or Ljt that it is a proton, kaon, or pion, respec- 
tively, are determined from the information provided by the 
hadron identification system. A track is identified as a proton 
if Lp/{Lp + Lk) > 0.6 and Lp/ {Lp + L„) > 0.6, as a kaon if 
Lk/ (Lk+Lti) > 0.6, or as a pion if L,i/(L/(- +L„) > 0.6. The 
efficiency for identifying a kaon (pion) is 85—95% depending 
on the momentum of the track, while the probability for a pion 
(kaon) to be misidentified as a kaon (pion) is 10—20%. The 
proton identification efficiency is 84% while the probability 
for a kaon or a pion to be misidentified as a proton is less than 
10%. 

We reconstruct A's from their decays to pK^. Each A can- 
didate must have a displaced vertex and the direction of its 
momentum vector must be consistent with an origin at the IP. 
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The proton-like daughter is required to satisfy the proton cri- 
teria described above, and no further selections are applied to 
the daughter tracks. The reconstructed A mass is required to 
be in the range 1.111 GeV/c^ < M^^- < 1.121 GeV/c^ H. 

Candidate mesons are reconstructed in the following 
two sub-decay channels: — > K^K^ and K^k^t\P, 
Tt" — > yy. The y's that constitute 7i'' candidates are required to 
have energies greater than 50 MeV if the y is reconstructed 
from the barrel ECL and greater than 100 MeV for the endcap 
ECL, and not be associated with any charged tracks in CDC. 

The energy asymmetry of y's from a ^-^j-^^^, is required 

to be less than 0.9. The mass of a Ti" candidate is required to be 
within the range 0.1 18 GeV/c^ < Myy < 0.150 GeV/c^ before 
a mass-constrained fit is applied to improve the T\P momentum 
resolution. We impose a cut on the invariant masses of the 
D" candidates, |M^-„+ - 1.865 GeV/c^j < 0.01 GeV/c^ and 
1.837 GeV/c2 < M^-„+„o < 1.885 GeV/c^ for D° K-%+ 
and — > K^%^%^, respectively, which retains about 87% of 
the signal. 

We reconstruct D*" mesons in the decay mode D*'-^ D^n^ 
with D° K^n^ only. Since the coming from the D*° 
decay is expected to have low energy, we adjust the photon 
selection criteria accordingly. The energy of y's that constitute 
Tt" candidates from a D*^ must be greater than 50 MeV. The 
energy asymmetry of the two y's is required to be less than 
0.6 and the di-photon invariant mass should be in the range 
0.120 GeV/c2 < Myy < 0.158 GeV/c^. For the D*'' candidates, 
we require 0.139 GeV/c^ < AM < 0.145 GeV/c^, where AM 
denotes the mass difference between D*" and D". 

Candidate B mesons are identified with two kinematic vari- 
ables calculated in the center-of-mass (CM) frame: the beam- 
energy-constrained mass Mbc = y ^beam ^ 

pI, and the energy 

difference AE = Eg — fibeam. where fibeam is the beam energy, 
and pb and Eg are the momentum and energy, respectively, of 
the reconstructed B meson. In order to reduce the contribu- 
tion from combinatoric backgrounds, we define the candidate 
region for B pAD" (pAZ)*") as 5.2 GeV/c^ < Mbc < 5.3 
GeV/c2, -0.1 GeV < A£ < 0.4 GeV and M^a < 3.4 (3.3) 
GeV/c^, where Mp\ denotes the invariant mass of the baryon 
pair. The lower bound in AE is chosen to exclude backgrounds 
from multibody baryonic B decays. From Monte Carlo (MC) 
simulations based on GEANT llSll . we define the signal re- 
gion as 5.27 GeV/c2 < Mbc < 5.29 GeV/c^ and |A£| < 0.05 
GeV. 

The dominant background for B 'pAlf' in the candi- 
date region is from continuum e+e^ qq {q ^ u, d, s, c) 
processes. We suppress the jet-like continuum background 
relative to the more spherical BB signal using a Fisher dis- 
criminant that combines seven event-shape variables derived 
from modified Fox- Wolfram moments |19i] as described in 
Ref. ll20ll . The Fisher discriminant is a linear combination of 
several variables with coefficients that are optimized to sep- 
arate signal and background. In addition to the Fisher dis- 
criminant, two variables cos 9^ and Az are used to form signal 
and background probability density functions (PDFs). The 
variable Qb is the angle between the reconstructed B direc- 



tion and the beam axis in the CM frame, and Az is the dif- 
ference between the z positions of the candidate B vertex and 
the vertex of the rest of the final state particles, presumably, 
from the other B in the T(4S) decay. The products of the 
above PDFs, obtained from signal and continuum MC simu- 
lations, give the event-by-event signal and background like- 
lihoods, Ls and Lg. We apply a selection on the likelihood 
ratio, 2^ = Ls/{Ls + Lb) to suppress background. Informa- 
tion associated with the accompanying B meson can also be 
used to distinguish B events from continuum events. The 
variables used are "q" and "r" from a B flavor-tagging al- 
gorithm ll2lll . The value of the preferred flavor q equals +1 
for B^/B^ and —1 for B^/B^. The B tagging quality fac- 
tor r ranges from zero for no flavor information to unity for 
unambiguous flavor assignment. Sets of ^ x r-dependent 'J{_ 
selection requirements are optimized by maximizing a fig- 
ure of merit defined as Ns / \/Ns+Nb, where Ns denotes the 
expected number of signal events based on MC simulation 
and the predicted branching fraction, and Nb denotes the 
expected number of background events from the continuum 
MC. The requirements on 2^ remove 75% (89%) of contin- 
uum background while retaining 88% (69%) of the signal for 
B pAD^ with D° K-%+ (D" R-k+k^). The contin- 
uum background suppression is not applied for B pAD*" 
since the optimal 3^ requirement is close to zero. 

In order to avoid multiple counting, in cases where more 
than one B candidate is found in a single event, we choose 
the one with the smallest = xl + Xa(+X^o)' where is 
calculated from the vertex fit to the B using 'p and mea- 
surements, the vertex fit to A using p and tracks, and the 
mass-constrained T\P fit if applicable. The fraction of events 
with multiple candidates are 2.3% (14.1%) of the sample ac- 
cording to MC simulations for B pAD" with — > K^%^ 
(D" ^ K-n+n\ and 17.8% for B pAD*^. The dominant 
background for B — > pAD*" is from Z?" pAD*^ cross- 
feed and B pAD*^ self cross-feed (both referred to as CF) 
events according to a MC simulation based on Pythia ||22|i . 
In CF events, two low-energy y's can form a n'^ candidate that 
is combined with a correctly reconstructed D", p and A from a 
B decay to form a candidate event in the signal region. These 
backgrounds cannot be distinguished from the signal in the 
Mbc — A£ two-dimensional fit alone, although their distribu- 
tions in Mbc and AE have a slightly wider spread than the sig- 
nal. We can, however, estimate this background contribution 
by analyzing the AM distribution in which signal events have a 
Gaussian shape and background events have a threshold func- 
tion shape as shown in Fig.[T] 

The signal yields the of B pAD'*''' modes are extracted 
from a two-dimensional extended unbinned maximum likeli- 
hood fit with the hkeUhood defined as 

!=1 J 

where is the total number of candidate events, Nj denotes 
the number of corresponding category events and P/ repre- 
sents the corresponding two-dimensional PDF in Mbc and AE; 
i denotes the i-th event, and j indicates the index of different 
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FIG 1: The AM distribution of the pAD*" MC sample with 

fit curves overlaid, where AM denotes the mass difference between 
D*° and Z)°. The solid curve is the overall fit result, the dashed curve 
shows the CF background and the black filled squares are the MC 
events. 

event categories in the fit. Thus, j could either indicate sig- 
nal or combinatorial background for the /JAZ)" case and in- 
cludes one more category (CF) for the pAD*" case. We use 
a Gaussian function to represent the signal and a dou- 
ble Gaussian function for the AE signal with parameters de- 
termined using MC simulations. Combinatorial background 
is described by an ARGUS function ll23ll and a second-order 
Chebyshev polynomial in the Mbc and AE distributions, re- 
spectively. 

Since it is difficult to separate the B pAD*'' signal and 
CF events in the fit, we estimate the number of CF events in 
the AM signal region (0.139 GeV/c^ < AM < 0.145 GeV/c^) 
from the fitted CF yield in the AM sideband region (0.15 
GeV/c2 < AM < 0.17 GeV/c^). The ratio of the area of the 
CF in the AM signal region to that in the sideband region 
is 26.0 ±0.9%, which is determined from MC samples of 
B pAD*^ (Fig. 1) and fiO ^ pAD*+ . The PDF used for 
CF events is a product of a Gaussian-like smoothed histogram 
for Mbc and a double Gaussian function for AE with param- 
eters determined using MC simulations. We fix the number 
of CF events in the Mbc ^ fit to determine the signal yield 
within the AM signal region. 

Figure |2] shows the result of the fit for B pAZ)". 
The fitted signal yields in the data sample are 26.5^5 g and 
35.6^[Q y events with statistical significances of 7.6 and 3.6 
standard deviations (a) for B — > pAD^, — > K^n+ and 
— > K^K^TiP, re spectively. The significance is defined as 
■\/ — 21n (Lo / Lynax ) , where Lq and Lmax are the likelihood val- 
ues returned by the fit with a signal yield fixed to zero and the 
nominal fit, respectively. The branching fractions are calcu- 
lated using the formula 

^ _ Nsignal 

where Nsig„ai, Ngg, e, and / are the fitted number of signal 
events, the number of BB pairs, the reconstruction efficiency, 
and the relevant sub-decay branching fractions; ® (A ^ 
pn-) ^ 63.9 ±0.5%, 'B{D^ ^ K-n+) =3.89±0.05%, 
'B{D'^ ^ = 13.9 ±0.5%, and ^dV) = 



TABLE I: Summary of the results: event yield, significance, effi- 
ciency, and branching fraction. 
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1.43+j];|±0.18 


B- 75AO*" 




2.1 


2.8 


1.53+',;^^ ±0.47 



6 1.9 ±2. 9% m. We assume that charged and neutral pairs 
are equally produced at the T(45'). 

To investigate the threshold enhancement feature, we de- 
termine the differential branching fractions in bins of M^a; 
the results obtained from the weighted averages of the fits to 
B pAD^\ D° K-%+ and -> R-Ti+Ti^ separately are 
shown in Fig. |3] where an enhancement near threshold is ev- 
ident. We fit the pA mass spectrum with a threshold func- 
tion and then reweight MC events to match the fitted thresh- 
old function in order to obtain a proper estimate of the recon- 
struction efficiency for signal events. The observed branching 
fractions are {l.39tlf^±Q.l6) x 10"^ forB" -^pAD^ with 
^K-%+ and (1.54;[J;^°±0.26) x 10"^ for B" -^pAD^ 
with K^ti^t6\ The weighted average of the branching 

fractions is (1.43+jJ-^^ ± 0.18) x 10"^ with a significance of 
8.1a, where the systematic uncertainties (described below) on 
the signal yield are also included in the significance evalua- 
tion. 

The fit results for B pAD*^ in the AM sideband region 
are shown in Fig. |4] (a, b). The number of CF events in 
the sideband is 11.6 ±5.4, which is used to estimate the 
number of CF events in the AM signal region, 3.0 ± 1.4, 
after scaling by the area ratio of CF (26. 0± 0.9%). We then 
fix the normalization of the CF component in the fit to the 
AM signal region [fit results are shown in Fig.|4](c, d)], and 
obtain a signal yield of 4.3^2 4 with a statistical significance 
of 2.2a. Assuming B — > pAD*^ and B pAD^ have 
the same pA spectrum, we determine «(B" pAD*^) 
to be (1.53;i;J^±0.47) x 10"^ In the absence of a sta- 
tistically compelling signal yield, we set an upper limit 
® (B- pAD*^) < 4.8 X 10"^ at the 90% confidence level 
using the Feldman-Cousins method ll24[l25ll . The information 
used to obtain the upper limit includes the number of events in 
the signal region (13) and 8.1 ± 1.4 background events. Here, 
the background that is integrated in the signal region, consists 
of 5.3±0.5 continuum events and 2.9±1.3 CF events from 
the fit to the AM sideband. The 11.7% additive systematic 
uncertainty due to the selection criteria is included in the 
determination of the upper limit on the branching fraction. 

Systematic uncertainties are estimated using high-statistics 
control samples. A track reconstruction efficiency uncertainty 
of 1.2% is assigned for each track. For the proton identifica- 
tion efficiency uncertainty, we use a A pK^ sample, and for 
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FIG 2; Distributions of A£ (a, c) for Mbc > 5.27 GeV/c^ and of Mbc 
(b, d) for \\E\ < 0.05 GeV; the top row is the fit result for B 
pAD", -s> K^TZ+ (a, b) and the bottom row for B ^> pAD*', 
— 7> K^Tf^Tt^ (c, d). The points with error bars are data; the solid 
curve shows the fit; the dashed curve represents the signal, and the 
dotted curve indicates continuum background. 




MpA (GeV/c^) 



FIG 3: Differential branching fraction (d® /dM-p/^) as a function of 
the pA mass for B pAD^*. Note that the last bin with the central 
value of 3 GeV/c^ has a bin width of 0.8 GeV/c^. The solid curve is 
a fit with a threshold function. 



K—K identification uncertainty we use a sample of kinemati- 
cally identified D*+ D^%+ , K'%+ decays. The aver- 

age efficiency discrepancy due to hadron identification differ- 
ences between data and MC simulations has been corrected 
for the final branching fraction measurements. The correc- 
tions due to the hadron identification are 10.7% and 10.6% for 
B pAD^ and B 'pAD*^, respectively. The uncertain- 
ties associated with the hadron identification corrections are 
4.2% for two protons (one from A decay), 0.5% for a charged 
pion, and 1.0% for a charged kaon. 

The Tt" selection uncertainty is found to be 5.0% by com- 
paring the ratios of efficiencies between — > K^K^ and 
Z)" — > K^Ti^n^ for data and MC samples. In the A recon- 
struction, we find an uncertainty of 4.1% from the differences 
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FIG 4: Distributions of A£ (a, c) for Mbc > 5.27 GeV/c- and of 
Mbc (b,d) for lA^I < 0.05 GeV; the top row is the fit result for B^ 
pAD*" in the AM sideband region (a,b) and the bottom row for i?^ — > 
pAD*" in the AM signal region (c, d). The points with error bars are 
data; the solid curve shows the result of the fit; the dot-dashed and 
dotted curve indicates the CF and continuum background; the dashed 
curve represents the signal. 



between data and MC for the efficiencies of tracks displaced 
from the interaction point, the A proper time distributions, 
and the A mass spectrum. The uncertainty due to the 'J{_ se- 
lection for B pAZ)", D'^ — > K^K^ is estimated from the 
control sample B D'^n^, /r^7t+7t^7t+ and is deter- 

mined to be 1.3%. The related uncertainty for B 'pAD^\ 

R-ii+TZ^ is 3.0% estimated from B 0°%-, D° 
A'^Tt+Tt". The uncertainties due to the £)" mass selection for 



Z)^ K-n+ and R-n+n^ are 1.9% and 1.6%, respec- 

tively. 

The dominant systematic uncertainty for B 'pAD^ is 
due to the modeling of PDFs, estimated by including a B ^■ 
pAD'^Tt or a nonresonant B — > 'pAK^K^ {K^h^t6^) compo- 
nent in the fit, modifying the efficiency after changing the sig- 
nal MpA distribution, and varying the parameters of the signal 
and background PDFs by one standard deviation using MC 
samples. The modeling uncertainties are 7.5% and 12.9% for 
B pAD^ with D" K-K+ and R-K+n^, respec- 

tively. The overall modeling uncertainty for B 'pAD*^ of 
28.6% is obtained from two kinds of PDF modifications. The 
parameters of the fixed CF component are varied by their ± lo 
statistical uncertainties, which were obtained from the fit to 
the AM sideband region. We also include an additional PDF 
for the combinatorial background based on the Pythia ||22|i 
b quark fragmentation process, e.g., B pAD", — >^ 



pA++D*'^, B 



■ pA^D*^, 



B 



plPD*'^, etc. 



The systematic uncertainties from the sub-decay branching 
fractions are calculated from the corresponding branching un- 
certainties in lH; they are 1.5% (3.7%) and 6.0% for B 



pAD", K-%+ (D^ K-n+n^) and B pAD 

spectively. The uncertainty in the number of BB pairs is 1 .4%. 
The total systematic uncertainties are 11.6% (17.1%) and 
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30.9% for B- pAD^ with D° K-n+ {D° R-n+n^) 
and B -> pAD*^, respectively. The final results are listed 
in Table m where the significance values are modified and in- 
clude the systematic uncertainty related to PDF modeling. 

In summary, using a sample of 657 x 10^ BB events, we 
report the first observation of B pAD^ with a branching 
fraction of (1.43+J5-^^ ± 0.18) x 10"^ and a significance 
of 8.1a. No significant signal is found for B 'pAD*^ 
and the corresponding upper limit is 4.8 x 10^^ at the 90% 
confidence level. We also observe a pA enhancement near 
threshold for B — > pAD", which is similar to a common 
feature found in charmless three-body baryonic B decays 
The measured B — > pAZ)" branching fraction agrees with 
the theoretical prediction of (1.14±0.26) x 10"^ lH. This 
indicates that the generalized factorization approach with pa- 
rameters determined from experimental data gives reasonable 
estimates for — > c decays. This information can be helpful 
for future theoretical studies of the angular distribution 
puzzle in the penguin-dominated processes, B — > p'pK^ 
and pAn^ ||5|l. The measured branching fraction for 

B —J- 'pAD^ can also be used to tune the parameters in the 
event generator, e.g., Pythia, for fragmentation processes 
involving b quarks. Although the current statistics for 
B pAZ)" are still too low to perform an angular analysis 
of the baryon-antibaryon system, the proposed super flavor 
factories ll26[|27ll offer promising venues for such studies. 
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